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INTRODUCTION 


New  spacecraft  designs  feature  larga  atructures 
characterised  by  low  natural  fraquanclaa  and  atrlngant  pointing 
and  vibration  requirements.  Thaaa  larga  apaca  structures  pose 
unique  and  difficult  control  problems.  These  problems  include 
system  bandwidth*  greater  than  structural  natural  frequencies; 
lack  of  accurate  information  about  the  dynamic  characteristic  of 
the  structure  being  controlled;  complicated  high-order  dynamics, 
including  non-linear  behavior;  and  stringent  requirements  for 
distributed  shape  control1.  An  important  part  of  the  solution  to 
these  control  problems  is  the  development  of  actuators  capable  of 
applying  force  or  torque  to  the  structures,  conventionally  these 
actuators  have  been  reaction  mass  actuators  or  distributed 
piezoelectric  materials. 

The  motivation  of  this  research  is  to  investigate  other 
innovative  actuator  designs.  In  particular,  the  direct  use  of 
electromagnetic  forces  is  seen  as  a  promising  approach. 
Originally  the  focus  of  this  research  pi'ogram  was  to  develop 
actuator  concepts  based  on  the  Lorentz  force  interactions  of 
current  carrying  wires.  These  concepts  range  from  discrete, 
distributed  coils  interacting  with  Earth's  magnetic  field  to 
distributed  winding  patterns  embedded  in  the  space  structure. 
These  actuators  are  expected  to  have  r  number  of  advantages 
including  high  bandwidth,  compatibility  with  composite 
structures,  and  low  mass.  In  addition  to  the  original  focus, 
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this  Phsss  I  program  expanded  its  foous  to  invsstigsts  ths 
applicability  of  aagnatostrictivs  matsrials  to  shapa  and 
vibration  control  actuators.  This  sxpandsd  focus  was  tha  raault 
of  tha  racant  davalopmant  of  high  parformanca  magnatostrictiva 
matarial  Tarfanol-D.  This  new  aatarial  offars  a  factor  of  10 
improvaaant  in  both  strain  capability  and  anargy  danaity  coaparad 
to  piezoelectric  matsrials . 

Tha  rasaarch  was  dividad  into  six  tasks.  First  a  flaxibla 
structura  modal  appropriata  for  control  and  actuator  assassmant 
was  devalopad.  Tha  naxt  two  tasks  invastigatad  actuator 
placamant  issuas  and  davalopad  dasign  dafinitions  for  tha 
advancad  actuator.  Tha  fourth  task  davalopad  controller 
algorithms.  The  fifth  task  avaluatad  the  parformanca  of  tha 
controlled  structura  using  tha  magnetic  coil  actuators.  Tha 
sixth  task  investigated  tha  affect  of  these  actuators  on  system 
issuas  such  as  interference  with  the  Barth's  magnetic  field. 

This  report  is  organized  as  follows.  The  first  major 
section  explains  the  model  development.  A  flexible  beam  model  of 
the  MAST  deployable  truss  space  structure  was  used  for  this 
study.  In  the  next  major  section,  which  is  the  main  emphasis  of 
this  research,  issues  pertaining  to  actuator  design  are  discussed 
and  an  initial  analysis  is  performed  on  each  design  method.  The 
magnetic  analysis  methods  needed  are  reviewed  briefly.  Five 
actuator  designs  are  identified.  The  first  is  the  interaction  of 
magnetic  coils  with  the  Earth's  magnetic  field.  The  second  is 
the  interaction  of  the  magnetic  coils  with  a  large  flux  source  on 
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another  part  of  the  spacecraft.  These  two  types  of  actuators  are 
shown  to  be  undesirable  for  the  MAST  application  yet  could  have 
advantages  for  control  of  other  types  of  flexible  structures. 
The  third  ie  the  interaction  of  distributed  Magnetic  coils  on  the 
flexible  bess  structure.  The  fourth  is  the  local  interaction  of 
magnetic  coils  placed  very  close  together.  The  final  approach  is 
the  use  of  a  new  magnetostrictive  material  called  Terfenol-D  to 
control  the  bean  through  strain.  Magnetostriction  is  the 
property  of  sone  materials  to  strain  under  the  presence  of  a 
magnetic  field.  The  design  and  analysis  of  these  approaches 
involved  investigations  into  system  issues  and  actuator  placement 
issues. 

In  the  next  major  section,  two  of  the  more  promising  methods 
the  distributed  coil  interaction  and  the  local  coil  interaction 
are  compared  in  more  detail.  Specifically,  the  control 
effectiveness  of  the  actuators,  which  is  based  on  excitation  mode 
and  actuator  location,  is  derived  and  compared.  The  local 
interaction  method  is  shown  to  be  more  effective  than  the 
distributed  on  both  a  per  mass  and  a  per  power  basis.  The  use  of 
Terfenol-D  is  also  a  very  promising  method,  as  is  shown;  however, 
more  comparison  analysis  needs  to  be  done.  It  is  proposed  that 
comparisons  of  the  local  interaction  and  the  magnetostrictive 
actuator  approaches  be  performed  as  an  initial  task  for  Phase  II 
research.  The  purpose  of  the  next  major  section,  Controller 
Design,  is  to  illustrate  how  the  local  interaction  actuators  may 
be  used  in  controlling  a  beam  structure  and  to  verify  the  control 
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effectiveness  concepts  that  vara  developed.  The  final  section  is 
the  conclusions  section. 


MODEL  DEVELOPMENT 

This  section  outlines  the  method  that  was  used  to  obtain  a 
large  flexible  structure  model  which  is  appropriate  to  the 
investigation  of  distributed  aagnetic  actuators.  The  model  used 
is  a  mathematical  representation  of  the  MAST  deployable  space 
truss  structure,  which  has  a  lowest  natural  frequency  of  0.2  Hz. 
This  section  contains  an  overview  of  the  model  requirements2,  a 
derivation  of  the  finite  element  equations  for  the  structure,  and 
finally  a  presentation  of  the  structure's  dynamic 
characteristics . 

Overview  of  the  Model  Requirements 

The  main  purpose  of  this  study  is  to  assess  the  advantages 
that  can  be  gained  by  using  distributed  magnetic  actuators  to 
achieve  structural  shape  and  vibration  control .  In  order  to  meet 
this  goal,  a  flexible  structure  model  is  needed  so  that  the 
closed-loop  performance  of  the  controlled  structure  with  the 
actuators  can  be  evaluated.  This  section  outlines  some  of  the 
desired  features  of  the  flexible  structu?.e  model. 

The  ideal  model  is  complex  enough  to  reflect  flexible 
structure  dynamical  features,  such  as  low  natural  frequencies  and 
lightly  damped  structural  modes.  On  the  other  hand,  the  model 
must  be  simple  enough  to  be  analytically  tractable  and  produce 
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useful  general  insights  into  system  behavior.  In  order  to  be 
compatible  with  the  control  software  tools  which  will  be  used, 
the  model  must  be  reducible  to  a  finite  state  format. 

To  satisfy  these  model  requirements,  a  two  degree-of-freedom 
finite  element  cantilevered  beam  model  subject  to  transverse 
shear  forces  and  bending  moments  was  chosen  for  this  study.  A 
beam  was  chosen  because  it  can  be  an  appropriate  simplified  model 
for  many  different  types  of  flexible  structures;  for  example,  a 
beam  could  be  a  simple  model  of  a  robotic  arm  or  a  deplt  yable 
truss  structure.  The  actual  parameters  for  this  beam  model  are 
compatible  with  the  MAST  structure.  This  structure  is  a  60  m 
long  deployable  truss  structure  which  was  designed  to  be  used  on 
space  shuttle  missions3.  Figure  1  is  an  artist's  conception  of 
MAST4.  Because  the  shuttle  is  much  more  massive  than  the 
structure,  a  cantilever  beam  is  a  good  model  of  the  system. 
Although  this  NASA  program  is  now  cancelled,  similar  flexible 
structures  are  likely  for  implementation  in  future  space 
programs.  The  finite  element  equations  of  the  beam  were 
manipulated  to  obtain  a  finite  state  representation  of  the  model. 

The  Finite  Element  Method  for  a  General  Structure 

The  finite  element  method  is  based  on  the  principle  of 
virtual  work  which  states  that  equilibrium  of  a  body  requires 
that  for  any  compatible,  small  virtual  displacements  (which 
satisfy  the  geometric  boundary  conditions) ,  the  total  internal 
virtual  work  is  equal  to  the  total  external  virtual  work. 
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Figure  1 
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Artist's  Conception  of  MAST  Structure 

Therefore,  for  a  general  three-dimensional  body  such  as  shown  in 
Figure  2,  the  principle  can  be  stated  as  follows: 

f  U'  fB  dV  +  f  US/  fs  dS  +  Si  Ui#  Fi 


t'r  dV  = 


(0 


where  the  overbar  denotes  virtual  quantities ,  and 

fB  -  [  fxB  fyB  f2B  ] '  =  Externally  Applied  Body  Forces 

fs  *  [  fxs  fys  f2s  ] #  =  Externally  Applied  Surface  Forces 

F*  ■  [  Fx*  Fy^  F2i  ]  •  =  Externally  Applied  Concentrated 

Forces 

U'  ■  [  U  V  W  ]'  «  Body  Displacements 

e  *  [  fxx  eyy  *  zz  7xy  7yz  'iz'x.}'  ~  Body  Strains 

f  **  t  f  xx  ^ yy  ^zz  ^xy  ^yz  f  zxl  *  B  Body  Stresses 

V  =  Volume 

S  *  Surface  Area 

Although  Equation  (1)  is  an  expression  of  equilibrium,  it  can 
also  satisfy  constitutive  and  compatibility  requirements  if  the 
problem  is  formulated  properly6. 

The  general  procedure  of  any  finite  element  problem  begins 
with  approximating  the  body  as  an  assemblage  of  discrete  finite 
elements  which  are  interconnected  at  nodal  points.  The  next  step 
is  to  obtain  individual  element  equations  which  satisfy  the 
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constitutive  and  compatibility  requirements7. 

Element  stresses  can  be  related  to  element  strains  by  the 
elasticity  matrix  (Equation  (2)).  When  this  relationship  is 
incorporated  into  the  principal  of  virtual  work,  constitutive 
requirements  are  satisfied8. 

rm  .  C®  +  fI®  (2) 

Cm  ■  Elasticity  Matrix  for  Element 
r1®  =  Element  Initial  Stress  Vector 
e®  ~  Element  Strain  Vector 

For  compatibility  requirements,  first  a  displacement 
interpolation  matrix,  H®,  is  formulated  to  relate  internal 
element  displacements  to  displacements  at  nodal  points  (Equation 
(3)).  Strain-displacement  relations  can  then  be  satisfied  by 
appropriate  manipulations  of  the  displacement  interpolation 
matrix  (Equation  (4))9. 

u®(x,y,z)  -  H® ( x , y , z )  u  (3) 

u®  =  Element  Displacement  Vector 
U  «  Vector  of  Global  Displacements  at  Nodal  Points 

em(X,y,z)  =  B® ( x , y , z )  U  (4) 

B®  =  Strain  Displacement  Matrix 


The  last  step  of  the  finite  element  method  is  to  assemble 
the  individual  element  equations  in  order  to  satisfy  the 
principle  of  virtual  work. 


Sm  «*'  dvn  =  2m 

U®'  fB®  dv®  +  Z 

U®s •  fSm  dsm  +  ui/pi 

J  v® 

J  v® 

Js® 

8 


(5) 


In  the  above  expression,  the  element  contributions  to  the  virtual 
work  are  calculated  by  integrating  over  the  elemental  volume,  V®, 
or  surface  area,  S®,  appropriately;  and  the  contributions  of  all 
the  elements  are  summed  up  to  obtain  the  total  virtual  work.  By 
substituting  Equations  (2),  (3),  and  (4)  into  Equation  (5), 
compatibility,  constitutive,  and  equilibrium  requirements  can  be 
satisfied10. 


U'  [Em  3®'C®  B®  dV®  ]  U 

J  V® 


U'  [ 


+ 


< 

(  Sm 

{  2m 


(H®'fB®  dVm  j 

ym 

1HSm#  fSm  dsm  j 
sm 

IB®'  r1  dV®}  +  F  ] 

V® 


(6) 


Since  the  principle  of  virtual  work  is  satisfied  for  any 
arbitrary  virtual  displacements,  which  do  not  violate  geometric 
boundary  conditions,  all  of  the  virtual  displacements  can  be  set 
equal  to  one.  Equation  (6)  can  now  be  rewritten  as: 

K  U  =  R  (7) 

K  -  B®'  C®  B®  dv®  “  Zjq  K® 

J  V® 

R  «  Rg  +  Rg  -  Rj  +  Rq 

RB  -  \  \  H®  ' f ®®  dV®  =  Zm  Rb® 

Jyffi 
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‘fffT  tom* HW»***i 


KtTfi.  ,  «r<>  W  ■  - 


Rs 

RI 

Rc 


-  m  f  Hs®' 

Js® 

-  m  f  B®' 

J  y® 


fs®  ds®  -  eb  rs® 


r1  dV® 


sn  Rj» 
2a  RC® 


D'Alembert's  principle  can  be  invoiced  to  incorporate  dynamic 
effects  into  the  model  (Equation  (8)  and  Equation  (9)). 


rB 


I. 


sm  I  H®'  (  fBm  “  Pm  H®  U  )  dV® 

I  V-H 


(8) 


U 


,m 


dU2  /  dt2 

Mass  density  of  element  m 


"■■•I, 


(9) 


p®  H®'  H®  dV®  =  E  M® 

I  yttl 

Therefore,  the  finite  element  equations  for  a  general  body  are11: 

M  U  +  K  U  «  R  (10) 

Boundary  conditions  are  imposed  by  splitting  up  the  problem 
into  known  displacements  and  unknown  displacements  and  then 
expressing  the  unknown  displacements  in  terms  of  the  known 
displacements  (Equations  (ll)  and  (12)).  The  resultant  set  of 
equations  has  the  same  ''spring-mass”  form  as  the  original  set  of 
equations,  only  the  load  vector  is  modified  in  such  a  way  that 
boundary  conditions  are  satisfied. 


M 


aa 


Mab 


Mba  Mbb 


Ub 


+ 

+ 


Kaa 

Kba 

10 


Kbb 


*1 1 

Ua 

- 

Ra 

_ub_ 

m 

Rb_ 

(ID 


Ua  »  Prescribed  Displacements 
Uj;,  -  Unknown  Displacements 

Mbb  ub  +  Kbb  ub  “  Rb  "Mba  ua  Kba  ua  (12) 

Once  the  unknown  displacements  have  been  solved  for,  the  unknown 
loads  can  be  obtained  as  is  shown  in  Equation  (13). 

Ra  -  Maa  Ua  +  Hab  ”b  +  *aa  Ua  +  Kab  Ub  (13) 

Finite  Element  Equations  for  a  2 -degree-of -freedom  Cantilever 
Beam 

The  flexible  structure  model  used  in  this  study  was  a 
cantilevered  beam  subject  to  shear  forces  and  bending  moments.  A 
computer  program  was  developed  to  give  the  finite  element 
equations  for  a  cantilevered  beam  with  variable  length  elements, 
1  (see  Figure  3) .  The  beam  used  has  a  constant  cross-sectional 
area  as  shown.  The  element  is  described  by  axial,  transverse, 
and  rotational  displacements.  Axial  displacements  are  later 
dropped  from  the  final  equations  because  the  bending  behavior  is 
of  more  interest  than  the  axial  behavior.  Allowing  the  program 
to  have  a  variable  number  of  variable  length  elements  gives  some 
degree  of  control  over  the  comp ’.exit}'  of  analysis. 

A  35,  140  state  element  mocV.i  was  used  for  natural  frequency  and 
mode  shape  computations.  Models  with  fewer  elements  were  used 
for  controller  design 

The  finite  element  equations  for  the  cantilever  beam 
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neglecting  damping  can  ba  written  aa  follows: 

M  U  +  K  U  »  R  (14) 

N+l 

K  »  Z  K® 

1 

N+l 

N  -  Z  M® 

1 

U  -  [  Ui  Wi  »2  02  •••  UN+1  WN+1  ®N+1  y 

m  (Nodal  Displacement  Vector) 

R  -  [  Pi  »!  Fj  M2  ...  FN+i  MN+1  ) ' 

*  (Concentrated  Nodal  Force  Vector) 
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l  2  l  H 


Kl~2  W'l  N  fj+  j 


Boundary  Conditions 

Wj  ■*  0 

6l  ■  0 

Bacause  tha  beam  is  cantilevered,  tha  boundary  conditions  at  tha 
wall  ara  zaro  transvarse  displacement  (w^  -  0)  and  zaro 
rotational  displacement  (e*  -  0) .  Tha  effect  of  thasa  boundary 
conditions ,  when  tha  problem  is  split  into  Known  and  unknown 
displacements  as  discussed  earlier,  is  to  reduce  tha  order  of  tha 
set  of  aquations  by  2.  Body,  surface,  and  initial  stresses  ara 
assumed  to  be  zero  in  this  mod  tl.  SatCon  actuators  will  provide 
concentrated  control  forces  to  the  beam.  Thus  tha  concentrated 
force  vector  R  is  dependent  on  the  number  and  nodal  location  of 
the  actuators. 

As  explained  previously,  tha  usual  method  for  defining 
stiffness  and  mass  matrices  involves  choosing  interpolations  to 
relate  displacements  within  the  element  tj  displacements  of  the 
nodes.  In  the  model  used  for  this  study,  the  element  stiffness 
and  mass  matrices  were  derived  using  the  "exact"  displacements 
obtained  by  solving  the  static  beam  Equations  (15)  and  (16) 12. 

Axial  Benavior: 

d  du 

—  (EA  —  )  -  0  (15) 

dx  dx 

x  =  distance  along  element  i  and 

E  =  Young's  Modulus 

A  -  bh 

u(x  «  0)  *  u(NODE  i) 

14 


u(x  -  1) 
u(x,y) 


'  u (NODE  i+1) 
u(x)  -  ye(x) 


Bending  Behavior: 


d2 

dx2 


0 


(16) 


v(X  -  0)  ■  V(NODE  i+1)  and  w(*  ■  1)  "  w(NODE  i-*-l) 

I  «  Moment  of  Inart ia  -  bh3/i2  (aaa  Figure  3) 

The  element  stiffness  and  mass  matrices  (neglecting  axial 
behavior)  are  shown  in  Figure  4  (The  structure  of  the  stiffness 
and  mass  matrices  which  incorporate  axial  behavior  is  such  that 
the  transverse  displacement,  w,  and  the  rotational  displacement, 
e,  do  not  depend  on  the  axial  displacement,  u) .  The  stiffness 
and  mass  matrices  for  the  whole  structure  are  obtained  by  adding 
together  the  appropriate  matrices  of  all  the  elements. 

To  transform  to  a  first  order  finite  state  set  of  equations 
needed  for  the  control  software,  a  new  state  vector  (which 
neglects  axial  behavior)  is  defined  and  Equation  (14)  is 
rewritten  as: 

U 

U  -  -M 


U  »  [w2  e2  ...  Wji+i  «h+1  w2  e2  •••  WN+1  ®N+l3T 

Equation  (17)  includes  the  effect  of  the  boundary  conditions  of 
aero  displacement  and  rotation  at  the  wall. 


15 


r 


C  3 

C  h 
C  — 
C  3 

C  ! 

C 

C  3 
C  h 

C 

C  2  1 

C 

C 

C  h 
C  -  - 
C 

C  1 

c 

C  3 

C  h 

C - 

C 

C  2  1 


h 

2  6  1 


h  h 

2  3  1 


(a)  Element  Stiffness  Matrix 
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(b)  Element  Mass  Matrix 
Figure  4 
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Dynamic  Chanotiriitioi  of  the  Bmi  Modal 

Tho  beam  nodal  uaad  In  thin  study  oould  ba  a  sinplifiad 
nodal  for  a  variaty  of  flaxibla  structures,  such  as  a  robotic  arm 
or  a  daployabla  truss  structure.  Tha  baan  nodal  chosen  is  a 
nathenatical  representation  of  tha  MAST  space  structure  (Figure 
1) .  A  baan  of  length  60  u>  (107  ft)  with  a  lowest  nstur?l 
frequency  of  1.15  rad/  sec  (0.184  Hi)  was  chosen  for  this 
analysis.  To  obtain  these  values ,  the  physical  parameters  used 
are  displayed  in  Table  1.  A  247  kg  mass  (544  Ibm)  with  a  moment 
of  inertia  of  20  kg  m2  (475  lbm  ft2)  was  situated  at  the  tip  of 
this  model13.  A  35  element  model  was  used  in  the  calculation  of 
natural  frequencies  (Table  2)  and  mode  shapes.  Figure  5  is  a 
plot  of  the  first  four  mode  shapes.  The  lowest  mode  in  this  plot 
has  a  large  transverse  deflection  at  the  tip.  As  the  excitation 
frequencies  get  higher,  more  energy  is  required  to  move  the  large 
tip  mass  and  the  tip  increasingly  behaves  like  a  pinned  end. 

Table  1.  Physical  Parameters  of  Beam  Model 
E  »  Young's  Modulus  »  5.2  x  107  N/m2  (7542  psi) 

v  -  Poisson's  Ratio  *  0 

h  -  Beam  Height  *»  1.7  m  (5.6  ft) 

b  -  Beam  Width  «  1.42  m  (4.66  ft) 

L  *  Beam  Length  «  60  m  (197  ft) 

p  *  Mass  Density  «  1.95  Kg/m3  (0.122  lbm/ft3) 
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NORMALIZED  BEAM  DISPLACEMENT 


Table  2.  Lowest  Four  Natural  Frequencies  of 
35  Element  Modal. 

0.194  Ha 

1.93  Ha 

5.7  Ha 

11.7  Ha 


POSITION  ALONG  LENGTH  OF  BEAM  (m) 


Figure  5  Maflft  ShflPfg 
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ACTUATOR  DESIGN 


This  section  is  a  discussion  of  the  actuator  design  concepts 
which  were  investigated.  Because  many  of  the  acutator  designs 
involve  the  production  of  forces  and  torques  by  interacting  cur¬ 
rents,  the  first  part  of  the  section  is  a  review  of  the  electro¬ 
magnetic  analysis  methods.  The  next  section  describes  the  five 
actuator  approaches  examined  and  discusses  some  of  the  advantages 
and  disadvantages  associated  with  each  particular  approach. 
Magnetic  Analysis  Methods 


Magnetic  Forces 


Circuits  carrying  electrical  current  produce  forces  on  each 
other.  For  the  general  case  depicted  in  Figure  6,  this  force  can 
be  obtained  by  the  magnetic  force  law 


Mo 

fab  ~  (  T“  ) *a  *b 
4ir 


1 

Sllb  *  (dla  x  El) 
r2 


/ 


(18) 


where  is  the  force  exerted  by  current  la  on  current  lb*  The 
line  integrals  are  evaluated  over  the  two  circuits,  and  the 
vectors  <Ua  and  dlb  are  incremental  elements  which  point  in  the 
direction  of  the  current  flow.  The  constant  nQ  -  4w  x  10  “7  N/A2 
(8.99  x  10“®  lb/A2)  is  the  permeability  of  free  space14.  The 
magnetic  force  can  be  written  another  way  so  that  it  is  more 
obvious  that  Newton's  third  law  holds. 

In  general  these  integrals  cannot  be  evaluated  analytically. 
They  are  performed  numerically.  However  by  rewriting  Equation 
(18)  ,  we  can  express  the  force  as  an  interaction  of  the 
current  Mb"  with  the  field  of  the  current  "a". 


Zab  =  (p  Silb  x  Ba  ( 19 ) 

where 

Ba  -  (  )  la(b  --  (dla  *  n)  (20) 

4k  Ja  r2 

This  is  the  Biot-Savart  Law  for  the  force  between  two  current 
loops15.  Ba  is  the  magnetic  induction  due  ro  circuit  a  at  the 
position  of  element  dlab  and  ls  measured  in  units  of  Tesla.  This 
expression  can  be  evaluated  for  simple  geometries.  The  flux  of 
the  magnetic  induction  through  a  surface  S  is  defined  by16 

*  -  |  fi  •  4A  (21) 

Another  method  can  be  used  to  calculate  the  force  exerted  on 
circuit  "a"  due  to  circuit  "bH.  This  method  is  based  on  energy 
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methods.  By  assuming  a  small  virtual  translation  of  one  coil  and 
then  using  the  conservation  of  energy,  the  magnetic  force  in  any 
direction  x  can  be  written  as 

Eab  -  I«  Ib  (  |;  )  (22) 

where  M  is  the  mutual  inductance  between  the  two  circuits  and  can 
be  related  to  the  magnetic  flux  linking  the  circuits17. 

*  ba  -  M  lb  (23) 

In  physical  situations,  these  interacting  circuits  which  produce 
force  on  one  another  are  made  up  of  coils  of  wire  and  can  thus  be 
referred  to  as  coils. 

fardrielfl  Approximation 

At  distances  which  are  large  compared  to  the  size  of  a  coil 
(i.e.,  at  least  10  coil  radii  away),  far  field  approximations  can 
be  made  for  the  magnetic  induction  and  the  interaction  force.  In 
this  case,  the  magnetic  field  can  be  described  using  spherical 
coordinates  r,e,  for  a  coil  located  at  the  origin  (Figure  7)18. 

Mo  2m 

Br  »*( - )  cos  e 

4*  r3 

Mo  » 

Be  »  ( - -  )  sin  e  (24) 

4*  r3 

B*  -  0 

where  a  *  ISjj  is  the  magnetic  dipole  moment  of  the  current  loop 
and  n  is  a  unit  vector  perpendicular  to  the  plane  of  the  coil. 
Because  field  falls  off  as  the  cube  of  the  distance  from  the 
origin  r3,  placing  a  second  interacting  coil  closer  to  the  origin 
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2 


Figure  7  Spherical  Coordinate  System 

will  produce  stronger  interaction  forces  and  torques.  The 
interaction  force  and  torque  due  to  coil  1  located  at  the  origin 
with  its  magnetic  moment  in  the  z  direction  on  coil  2  located  at 
(R,e,4)  (see  Figure  8)  can  be  written  as  below19 

E  =  12  7  (Bl  •  D2)  -  (»2  *  7)  Bl  (25) 

L  =  *2  x  Bl  (26) 

Notice  that  the  interaction  force  is  a  function  of  the  gradient 
of  the  magnetic  field  and  therefore  is  proportional  to  the 
inverse  of  the  fourth  power  of  the  separation  distance  (F  a 
1/r4).  The  torque,  on  the  other  hand,  is  proportional  to  the 
inverse  of  the  third  power  of  the  separation  distance  (r  a  1/r3) . 
Thus  the  strength  of  the  force  falls  off  more  quickly  with 
distance  than  does  the  strength  of  the  torque.  Because  we  are 
considering  actuation  schemes  which  involve  both  force  and  torque 
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capability  and  these  capabilities  are  sometimes  traded  off  in  the 
design  process,  this  fact  will  have  important  implications 
regarding  the  actuator  effectiveness.  When  the  far  field 
approximation  is  used,  the  coils  are  sometimes  referred  as 
dipoles. 

Example  •‘Standard"  Dipole 

Host  of  the  actuator  approaches  examined  will  use  control 
coils.  The  remainder  of  this  report  will  make  extensive  use  of  a 
"standard  size"  control  coil  for  scaling  analysis.  This 
"standard"  coil,  shown  in  Figure  9,  has  an  outer  radius  (rQ)  of  5 
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cm  (1.97  in),  an  Innar  radius  (ri)  of  2.9  on  (0.98  in),  a  height 
(h)  of  10  cn  (3.94  in)  and  a  currant  danaity  (J)  of  10*  A/m2  (645 
A/in2) .  This  currant  density  is  vary  conservative,  typical  of 
household  wiring.  For  reference  purposes,  the  volume  of  this 
coil  is  589  cm3  (35.9  in3),  and  its  mass  is  1.6  kg  (3.5  lbm)  if 
made  of  aluminum. 

h-r.Hnh- 

*  !  ' 

•  *  * 

rQ  -  0.05  a 
rj.  -  0.025  m 
h  -  0.1  n 
J  -  106  A/m2 

Figure  9  "Standard". .  Size  Coil 


The  magnetic  moment  of  the  "standard"  dipole  can  be 
calculated  using  Equation  (27) . 

■  *=0.5  f  1  x  J  dV  =  whJ(rQ3  -  ri3)  =  11.4  A  m2  (123  A  ft2)  (27) 


The  interaction  force  and  torque  between  two  of  these  "standard" 
dipoles  can  be  calculated  using  Equation  (25)  and  (26) .  For  a 
configuration  shown  in  Figure  10,  in  which  the  dipole  moments  are 
perpendicular  to  each  other,  and  a  coil  separation  distance  of  1 
m  (3.28  ft)  the  resulting  force  and  torque  on  coil  2  are 
3  Pol»lll»2l 


X  «  3.9  X  10"5  N(8.8  X  10“6  lb)* 


(28) 


*  r' 
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Figure  10  Configuration  for  Example  Calculations  of 

interaction  Force  and  Torque  on  Two  Standard 
Size  Coils 


PoIbiIIB2I  _  K 

x.  =  - - - - x  =  2.6  X  10"5  N  m  (1.9  x  10“5  ft  lb)Y  (29) 

2ir  r3 

Note  that  for  this  example  configuration,  the  force  and  torque 
are  perpendicular  both  to  each  other  and  to  the  axis  connecting 
the  coils.  Importantly,  a  system  containing  both  dipoles  will 
see  no  net  force  or  torque,  although  each  dipole  experiences  both 
torque  and  force.  The  equilibrium  of  a  system  containing  both 
dipoles  is  a  direct  consequence  of  Newton's  third  law.  The  force 
and  torque  on  coil  1  can  be  found  using  Newton's  third  law.  For 
example,  the  force  on  coil  1  is  equal  in  magnitude  and  opposite 
in  direction  to  the  force  on  coil  2  by  force  equilibrium. 


Maxwell's  Mutual  Inductance  Formula  for  Circular  Filaments 

When  the  dipole  approximation  can  not  be  used,  another 
method  must  be  used  to  obtain  the  interaction  forces. 

For  two  equal  radius  circular  filaments  coaxially  located  as 
shown  in  Figure  11,  the  mutual  inductance  can  be  expressed  in  a 
formula  given  by  Maxwell20: 
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(30) 


M  »  4  t  a  {  (  |  -  k)  F  -  |  E  ) 

where 

k  «  2  /  Jl  + 

a  -  radius  of  circular  filaments 

r  ■  distance  between  centers  of  filaments 

F  ■  Complete  elliptic  integral  of  the  first  kind 
with  respect  to  modulus  k 

E  -  Complete  elliptic  integral  of  the  second  kind  with  respect 
to  modulus  k 


T 


Figure  ii  Two  Ccaxi^l.  ssrufll-radiu?  circular  EUaments 

By  differentiating  this  expression  with  respect  to  r  and 
multiplying  by  the  products  of  the  currents,  as  in  Equation  (22) , 
an  expression  for  the  force  can  be  written.  This  differentiation 
can  be  performed  numerically  by  using  the  polynomial 
approximations  for  the  elliptic  integrals  found  in  Reference  21. 
It  turns  out  that  the  interaction  force  between  the  two  current 
loops  is  a  function  of  the  ratio  of  the  separation  distance  to 
the  coil  radius  (r/a) .  The  magnitude  of  the  force  normalized  by 
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the  product  of  currents  is  shown  in  Figure  12.  This  plot  shows 
that  at  near  field  (r/a  <l)  the  filament  force  is  inversely 
proportional  to  r/a.  For  far  field  (r/a  >10) ,  the  force  goes  as 
the  negative  fourth  power  of  r/a.  This  can  be  verified  using  the 
dipole  far  field  approximations  for  force  as  shown  earlier. 


Magnitude  of  Filament  Force  v$  r/a 


r/a  (cm/cm) 

Figure  12  Heat  Field.  Filament.  intejrflsUon  Forge 

ACTUATOR  APPROACHES 

This  section  describes  the  actuator  design  approaches 
examined  in  this  study.  Most  of  these  designs  involve  the 
interaction  of  a  magnetic  coil  with  a  magnetic  field  from  another 
source.  Hence  the  actuator  designs  are  based  on  the  analysis 
developed  in  the  previous  section.  The  first  actuator  design 
uses  a  number  of  distributed  coils  at  various  places  in  the 
flexible  structure  that  interact  with  the  Earth's  magnetic  field 
to  produce  control  torques.  The  second  design  approach  uses  a 
large  magnetic  flux  source  and  produce  torques  on  the  structure 
by  interaction  of  distributed  coils  and  this  flux  source.  The 
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third  design  approach  products  forces  and  torquas  on  tha 
structure  by  tha  interaction  of  distributed  coile.  The  fourth 
design  approach  uses  a  local  interaction  scheme  that  involves 
placing  coils  close  together  on  structural  truss  producing  force 
or  strain  that  is  transmitted  to  the  overall  beam  structure.  The 
final  approach  involves  the  use  of  a  magnetostrictive  material  in 
order  to  control  the  structure. 

Interaction  of  Distributed  Coils  with  Barth's  Magnetic  Field 

This  actuator  design  uses  magnetic  coils  distributed  along 
the  structure  interacting  with  the  Earth's  magnetic  field  and 
producing  control  torques  on  the  structure.  This  actuator 
concept  is  illustrated  schematically  in  Figure  13. 

The  Earth's  magnetic  field  can  be  modelled  as  that  of  a 
dipole.  The  magnitude  of  this  field  can  change  for  a  variety  of 
reasons.  Solar  wind,  solar  radiation  and  sunspots,  for  example, 
can  alter  the  field's  strength.  The  dipole  strength  is 
decreasing  by  0.05%  per  year22.  The  field  intensity  also  varies 
with  longitude,  latitude,  and  altitude.  The  variation  of  the 
maximum  value  of  the  Earth's  magnetic  field  with  altitude  is 
shown  in  Figure  14.  Notice  that  as  can  be  expected  from  the  far- 
field  equations  for  magnetic  strength  of  a  dipole  (Equation 
(24)),  the  strength  of  the  field  decreases  with  altitude. 
Altitudes  typical  of  low-Earth-orbit  applications  are  50  km. 
This  includes  space  shuttle  applications  such  as  MAST.  At  this 
altitude  the  magnetic  field  strength  is  on  the  order  cf  0.5  Gauss 
(  50  microtesla) .  This  field  strength  will  be  used  in 
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Figure  14  Earth's  MflgDfifclS. Eiftlfl  vs  Alfcitufla 
the  calculations  that  follow. 

Because  the  Earth's  magnetic  field  is  relatively  constant 
(i.e.,  the  gradient  is  very  small),  this  actuation  system  is  not 
good  at  producing  control  forces  on  the  structure.  This  is 
because  the  dipole  force  is  proportional  to  the  gradient  of  the 
magnetic  field  (Equation  25) .  The  control  torques  on  the 
structure,  on  the  other  har.d,  can  be  calculated  from  Equation 
(26) .  The  maximum  torque  that  this  technique  can  produce  was 
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calculated  by  ut>ing  a  ralativaly  large  control  dipole  and 
assuming  an  orientation  relative  to  the  Earth's  magnetic  field 
which  maximises  the  torque.  A  "standard"  site  control  dipole  of 
(Figure  9)  was  used  to  calculate  the  torque  available  at  an 
altitude  typical  for  low-Earth-orbit  space  applications  (such  as 
the  shuttle).  At  this  altitude,  the  earth's  magnetic  field  is  50 
microtasla  (0.5  Gauss).  The  resulting  torque  and  torque  per  unit 
of  actuator  mass  are  presented  below. 

r  -  5.7  x  10**4  rn  (4.2  x  10”4  ft  lb)  (31) 

r  Nm 

-  -  3.5  xlO”4  —  (1.2  x  10“4  ft  lb/lbm) 

mass  Kg 

These  torques  are  rather  small  for  controlling  the  beam 
model  used  in  this  study.  Another  disadvantage  of  this  actuator 
design  for  control  of  the  MAST  model  is  its  dependence  on  the 
orientation  of  the  actuator  to  the  Earth's  magnetic  field.  A 
control  system  which  incorporates  these  actuators  for  shape  and 
vibration  control  would  have  to  sense  and  control  positioning  of 
the  actuator  with  x;es>>ect  to  inertial  space. 

Interaction  of  Distributed  Coils  with  a  Large  Flux  source 

This  actuator  design  concept  is  similar  to  the  above  except 

that  instead  of  the  Earth's  magnetic  field,  a  large  flux  source 

¥ 

is  mounted  on  the^  spacecraft.  The  field  produced  by  this  flux 
source  interacts  with.,  distributed  coils  and  produces  control 
torques.  A  schematic  of  this  approach  is  shown  in  Figure  15. 

A  superconducting  solenoid  with  a  current  density  of  10s 
A/m2  (9.3  x  10®  A/ ft2)  magnetic  moment  of  107  Am2  (l  x  10®  A  ft2) 
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was  uaad  as  tha  baaalina  flux  aourca .  Tha  dimensions  of  this 
solanoid  ara  shown  in  Figura  16.  This  is  a  large  solanoid.  Its 
height  (h)  is  1  ■  (3.3  ft)  and  its  outar  radius  (r©)  is  0.5  m 
(1.6  ft).  Its  volume  is  0.6  a3  (21  ft3)  and  its  nass  is 
approximately  5200  kg  (1.1  x  104  lbs) . 


Figura  15  AstUfttor  Approach  2; _ intaraction  of  Distributed 

Colls  with  Laraa  Flux  Sourca 
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rQ  *  0.5  m 

ri  =  0.25  m 

h  -  1  a 

j  =  io8  A/a2 

Figure  16  Large  Flux  Source.  Solenoid  Dimensions 

At  a  distance  of  50  a  (164  ft),  the  maximum  torque  on  the 
"standard"  size  dipole  (Figure  9)  is  1.05  x  10~4  Nm  (8  x  10~5  ft 
lb) .  As  in  the  previous  actuator  approach,  this  torque  is 
dependent  on  the  orientation  of  the  control  dipole  with  respect 
to  the  large  flux  source.  The  disadvantage  of  this  approach  is 
that  torques  can  only  be  produced  in  two  directions.  This  is 
clear  from  the  dipole  equation  for  torque  (Equation  (26))  which 
states  that  the  torque  goes  as  the  cross  product  of  the  magnetic 
moment  and  the  magnetic  field;  thus,  torque  cannot  be  produced  in 
a  direction  which  is  parallel  to  the  magnetic  moment  of  the  large 
flux  source.  Another  disadvantage  for  application  of  these 
actuators  to  the  MAST  structure  is  that  the  large  flux  source 
experiences  a  large  torque  due  to  the  Earth's  magnetic  field. 
This  torque  is  about  570  Nm  (420  ft  lb)  at  an  altitude  place  when 
the  Earth's  magnetic  field  strength  is  50  microtesla  (0.5  Gauss) 
(i.e.,  an  altitude  typical  for  low-Earth-orbit  shuttle 
applications) .  Because  MAST  is  designed  to  be  deployed  from  the 
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shuttle  payload  bay,  this  torque  would  have  to  be  carried  by  the 
shuttle.  Flexible  structure  control  of  the  MAST  beam  would  have 
to  be  coupled  with  attitude  control  of  the  shuttle  In  order  for 
this  actuation  scheme  to  be  used.  Thus  this  actuator  approach  Is 
not  desirable  for  the  MAST  application.  This  scheme,  on  the 
other  hand,  may  be  more  advantageously  employed  for  another  type 
of  spacecraft  that  combines  attitude  control  with  flexible 
structure  control.  In  any  case,  the  large  flux  source  would  need 
to  be  superconducting  for  acceptable  power  consumption. 

Interaction  of  Distributed  Coils 

This  actuator  design  approach  places  a  number  of  magnetic 
coils  along  the  structure  and  obtains  forces  and  torques  by  the 
dipole  interaction  that  exists  between  them.  A  schematic  of  this 
approach  is  shown  in  Figure  17. 

Three  configurations  are  considered  which  encompass  the 
spectrum  of  all  possible  configurations  (i.e.,  a  general 
configuration  would  be  a  linear  combination  of  these  three) 
(Figure  18) .  The  first  dipole  is  placed  at  the  origin  with  its 
magnetic  moment  in  the  z  direction,  and  the  position  of  the 
second  dipole  is  varied  for  the  three  different  configurations. 
In  the  first  configuration,  the  second  dipole  is  offset  from  the 
first  by  a  distance,  r,  along  the  y  axis  and  has  a  magnetic 
moment  in  the  z  direction.  In  the  second  configuration,  the 
second  dipole  is  offset  a  distance  r  along  the  z  axis  and  its 
magnetic  moment  is  in  the  x  direction.  The  third  configuration 
has  the  second  dipole  on  the  z  axis  as  well,  but  with  the 
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magnetic  moment  in  the  z  direction.  The  dipole  forces  and 
torques  are  computed  as  a  function  of  magnetic  moment  using  the 
far-field  equations  developed  in  the  previous  section  (Equations 
(25)  and  (26))  and  Newton's  third  law.  The  configurations  and 
force/torque  on  coil  2  equations  are  shown  in  Figure  18.  The 
force  in  all  three  con-  figurations  is  proportional  to  1/r4  as 
expected  from  an  earlier  discussion  about  Equation  (25) .  The 
force  is  strongest  in  the  third  configuration  by  a  factor  of  2. 
This  configuration,  however,  has  no  torque  capability.  Only  the 
second  configuration  is  capable  of  producing  both  force  and 
torque,  as  can  be  seen  in  the 

35 


Configuration  2 
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Figure  18  (b)  Configuration  of  Dlstriouted  Coils  and  Tti 
Associated  Dipole  Interaction  Force7Torcnj 
Expressions  for  Force/Toroue  on  Coil  2 
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Configuration  3 
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equations  in  Figure  18.  The  magnitude  of  this  torque  is 
proportional  to  the  inverse  of  the  cube  of  the  separation 
distance  (r  a  1/r3) .  This  was  also  established  earlier  by 
consideration  of  the  dipole  Equation  (26) .  The  direction  of  the 
force  and  torque  capability  in  the  second  configuration  is 
appropriate  for  the  control  of  beam  bending  behavior  (Figure  19) . 
This  figure  shows  both  the  forces  and  the  torques  are  in 
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directions  which  support  control  of  bending  about  the  y  axis. 
The  sum  of  these  forces  and  the  sum  of  these  torques  is  zero  as 
can  be  expected  from  Newton's  Third  Law  of  Action-Reaction. 

Calculations  were  performed  to  determine  the  magnitude  of 
the  force  and  torque  capabilities  as  a  function  of  separation 
distance,  r,  and  current  density,  J.  Figures  20,  21,  and  22  are 
plots  of  coil  2  forces  and  torques  normalized  by  the  product  of 
the  magnetic  moments  of  the  two  coils  vs.  dipole  separation  dis¬ 
tance.  These  plots  indicate  the  effect  of  separation  distance. 
As  established  earlier  by  consideration  of  the  dipole  equations 
(Equations  (24),  (25),  (26)),  the  strength  of  the  force  or  torque 
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Figure  20  Force  vs  Separation  Distance  for  ConfiguratlP-D 

is  weaker  as  the  distance  between  the  dipoles  increases. 

Comparison  of  these  three  configurations  show  that  the  force 
is  strongest  when  the  dipole  moments  are  aligned  with  each  other 
(the  third  configuration).  If  the  dipoles  have  a  magnetic  dipole 
moment  of  11.4  A  ir.2  (123  A  ft2)  as  in  the  "standard"  dipole,  this 
force  capability  is  on  the  order  of  10~4  Newtons  (10~5  lb)  for  a 
aparation  distance  of  1/2  m  (1.6  ft).  The  other  two 
configurations  give  forces  that  are  half  as  strong  for  the  same 
magnetic  moments  and  separation  distance.  The  second 
configuration,  however,  is  the  only  one  which  is  capable  of 
producing  a  torque.  This  is  significant  as  the  torque  falls  off 
as  1/r3  while  the  force  falls  off  as  1/r4  as  shown  earlier. 

effect  of  current  density  is  shown  in  Figures  23,  24, 
and  25.  Plotted  is  coil  2  force  or  torque  vs  current  density. 
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"Standard"  size  dipoles  (Figure  (9))  are  used  with  a  separation 
distance  of  1  m  (3.3  ft)  (a  *  11.4  Am2  (123  A  ft2)).  Because  the 
magnetic  moment  is  proportional  to  the  current  density,  the 
torques  and  forces  are  proportional  to  the  square  of  the  current 
density.  Typical  values  of  current  density  range  from  106  A/m2 
(9  x  10*  A/ft2)  in  household  wiring  applications  *.o  10®  A/m2  (9  x 
106  A/ ft2)  in  superconductors.  The  plots  show  values  on  the  order 
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Figure  25  Force  vs  Current  Density  for 
Configuration  3 

of  0.4  N  (0.09  lb)  and  0.25  Nm  (0.18  ft  lb)  on  coil  2  for  the 
second  configuration  for  a  current  density  of  108  A/m2  (9  x  104 
A/ft2).  The  third  configuration  has  a  force  which  is  twice  as 
strong  (0.8  N)  (0.18  lb)  for  this  current  density  but  no  torque. 
Claims  of  current  densities  of  1010  A/m2  (9  x  108  A/ ft2)  have 
been  reported  for  new  high-temperature  superconductors,  raising 
the  force  and  torque  potential  of  this  actuator  system  by  a 
factor  of  10000. 

Further  investigation  of  the  feasibility  of  this  approach  is 
provided  in  a  following  section. 


local  Interaction  of  Magnetic  Coils  on  a  Trues  Section 

As  established  earlier,  stronger  forces  and  torques  can  be 
obtained  by  placing  the  coils  closer  together  (Equations  (25)  and 
(26)).  The  fourth  actuator  approach  is  different  from  the 
previous  concepts  in  that  the  coils  are  designed  to  be  very  close 
to  each  other  to  take  advantage  of  the  higher  force  capability. 

This  approach  is  to  place  magnetic  coils  on  a  truss  section 
of  the  MAST  beam  model.  A  schematic  of  this  approach  is  shown  in 
Figure  26.  Notice  that  these  coils  are  in  the  configuration 
where  their  magnetic  moments  are  aligned  on  the  same  axis.  The 
dipole  equations  (Equations  (25)  and  (26))  predict  a  force  along 
this  axis.  This  configuration  is  the  same  as  configuration  3  in 


BEAM  WITH  TRUSS  SECTION 


Figure  26  Astuatttr  ABBta&Sh  4; _ Local  Coil  Interaction 
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the  previous  discussion.  However,  unlike  the  previous  actuator 
approach,  the  dipole  equations  (Equations  (25)  and  (26))  can  not 
he  used  in  this  approach  because  they  break  down  at  close 
separation  distances  (r  <  10  coil  radii) .  Thus  the  alternate 
methods  using  Maxwell's  mutual  inductance  formula  (Equation  (30)) 
will  be  necessary  for  the  analysis. 

In  addition  to  more  complicated  methods  to  compute  the  coil 
force,  the  local  interaction  actuator  approach  requires  more 
complicated  modelling  for  the  MAST  beam-actuator  interaction.  As 
shown  in  Figure  27,  there  are  two  extreme  ways  to  model  this 
interaction  with  the  beam.  One  way  is  as  if  the  beam  is  much 
more  massive  than  the  truss  and  acts  essentially  as  a  kinematic 
constraint  preventing  met ion.  In  this  case,  the  force  at  the  end 
of  the  truss  is  completely  transmitted  to  the  beam.  If  the  coils 
are  centrally  located  in  the  truss,  and  the  truss  is  modelled  as 
a  continuous  spring,  the  force  transmitted  to  the  beam  is  coil 
force  reduced  by  a  factor  which  is  the  ratio  of  the  separation 
distance  (r)  to  the  truss  length  (1) .  Thus  there  is  an  inherent 
tradeoff  between  placing  the  coils  close  together  to  obtain  a 
stronger  magnetic  force  and  placing  them  close  to  the  end  of  the 
truss  so  that  more  force  is  transmitted  to  the  overall  beam 
structure . 

In  the  other  extreme  model,  the  beam  is  assumed  to  deflect 
due  to  the  magnetic  coil  interaction  force.  In  this  model,  no 
force  is  transmitted  to  the  beam.  As  emphasized  in  the  Figure 
27,  the  deflection  of  the  beam  under  this  modelling  approximation 
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Figure  27  Extreme  Models  of  Beam-Truss  Interaction  for 

tocal  Actuator  ABpraflsh 


is  also  proportional  to  the  ratio  of  the  separation  distance  (r) 
to  the  truss  length  (1) .  Thus  the  tradeoff  between  placing  the 
coils  close  together  and  separating  them  is  also  in  effect  here. 
The  two  modelling  extremes  can  be  used  to  find  two  extreme 
approximations  of  what  will  actually  happen.  In  actuality,  there 
will  be  some  force  and  some  motion  transmitted  to  the  beam. 

The  advantage  of  the  local  interaction  approach  is  that 
because  the  coils  are  closer  together,  they  are  able  to  produce  a 
stronger  force  than  the  other  type  actuator  designs.  To  analyze 
this  approach,  the  near  field  calculation  of  the  mutual 
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inductance  in  Maxwell's  formula  (Equation  (30))  is  necessary. 
This  expression  is  then  differentiated  to  find  the  interaction 
force  (Equation  (22)).  As  stated  earlier,  polynomial 
approximations  for  the  elliptic  integrals  in  the  expression  for 
mutual  inductance  and  numerical  differentiation  was  used  to 
calculate  the  coil  interaction  force.  The  results  were  shown  in 
Figure  12,  a  plot  of  force  normalized  by  the  product  of  the  coil 
currents  vs  the  ratio  of  coil  separation  to  coil  radius  (r/a) . 
This  plot  can  be  used  to  find  the  magnitude  of  the  force  for 

given  values  of  separation  distance,  coil  radius,  and  coil 

1 

currents . 

An  example  calculation  that  is  applicable  to  the  MAST  beam 

is  included  to  illustrate  this  approach.  For  reference,  some 

MAST  trusses  have  a  diameter  of  approximately  2.54  cm  (1  inch). 

This  value  is  used  for  the  diameter  of  two  coils  centered  in  such 

a  truss  in  a  configuration  so  that  their  magnetic  moments  are 

aligned  with  the  axis  of  the  truss  as  in  Figure  26.  The  coils 

are  assumed  to  be  separated  by  a  distance  of  1.27  cm  (1/2  inch) 

so  that  r/a  is  equal  to  1.  A  current  density  108  A/m2  and  a  wire 

cross-sectional  area  of  0.25  cm2  (3.9  x  10”2  in2)  is  assumed  for 

» 

both  coils.  Under  these  conditions  the  coils  produce  an 
interaction  force  on  each  other  of  approximately  8  x  10“7  x 
(108) 2  x  (2.5  x  10-5)2  =  5  n  (1.1  lb).  This  is  over  a  factor  of 
5  increase  in  comparable  calculations  for  the  distributed  coil 
interaction  case.  Because  both  the  distributed  and  the  local 
coil  interaction  actuator  designs  can  be  used  in  the  MAST 
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structure,  these  will  be  compared  in  further  detail  in  the  next 
major  section  of  this  report. 


SatCon  has  expanded  the  focus  of  this  research  program  to 
include  an  actuator  approach  which  uses  a  new  magnetostrictive 
material,  Terfenol-D.  This  approach  shows  much  potential  as  an 
actuator  for  intelligent  structures.  While  the  other  approaches 
are  based  on  magnetic  forces  which  arise  due  to  the  interaction 
of  currents,  this  approach  is  based  on  the  property  of 
magnetostrictivity.  This  means  the  material  will  produce  strains 
under  the  influence  of  a  magnetic  field.  This  actuator  approach 
is  now  very  promising  because  the  new  material,  Terfenol-D,  has 
shown  much  stronger  magnetostrictive  properties  than  other 
materials.  Furthermore,  it  has  a  factor  of  10  improvement  in 
both  strain  capability  and  energy  density  of  piezoelectric 
materials.  SatCon  has  included  some  preliminary  results  in  this 
report.  Further  analysis  is  proposed  under  an  initial  task  in 
Phase  II,  to  compare  the  use  of  Terfenol-D  with  the  other 
actuator  approaches  in  this  program.  Phase  II  will  then  identify 
one  concept  to  develop  into  a  baseline  prototype  lab  model. 

Magnetostriction 

Useful  magnetostrictive  materials  produce  large  strains 
(extensions)  when  subjected  to  magnetic  fields.  in  most 
materials  this  effect  is  quite  small.  For  example, 
magnetostrictive  strains  in  Nickel  are  on  the  order  of  10  parts 
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per  million  (microetrain) .  The  rare  earth  elements,  however, 
possess  a  number  of  extraordinary  magnetic  properties.  One  of 
these  properties  is  the  ability  to  produce  large  magnetostrictive 
strains.  During  the  1960 's  and  early  70 #s  a  variety  of  rare 
earth  materials  with  large  magnetostrictive  strains,  on  the  order 
of  1%  (10,000  microstrain)  ware  found.  These  materials,  however, 
only  exhibited  these  large  magnetostrictive  strains  at  low 
temperatures . 

During  the  1970 's,  a  variety  of  new  rare  earth  compounds 
that  exhibit  high  levels  of  magnetostriction  at  room  temperature 
were  developed.  These  rare  earth/ iron  compounds  have  large 
^oagnetostrictive  strains  of  up  to  0.2%  (2000  microstrain)  and 
Curie  temperatures  of  over  500  degrees  Fahrenheit.  Researchers 
at  the  Ames  Laboratory  of  the  Department  of  Energy  and  at  the 
Naval  Surface  Weapons  Center  have  further  developed  these 
compounds  and  their  manufacturing  processes  to  the  point  where 
they  have  recently  become  commercially  available  under  the  brand- 
name  Terfenol-D.  Terfenol-D  is  formed  from  the  rare  earths 
terbium  and  dysprosium  and  from  iron  (Tb^Dy^Fe.^  g3)  which  is 
directionally  solidified  by  free-stand-zone-melt  (FSZM)  process 
to  near  single  crystal  or  directionally  solidified  by  the 
Modified  Bridgman  process.  The  resulting  material  is  available 
in  rods  of  up  to  1.5  inches  (3.8  cm)  in  diameter.  This  stock  can 
then  be  sliced  into  lamination  thicknesses  if  desired. 

Because  of  these  recent  developments,  large-strain 
magnetostrictive  materials  capable  of  room  temperature  operation 


50 


are  now  available  as  the  engineering  material  Terfenol-D.  The 
most  important  property  of  Terfenol-D  is  its  magnetostrictive 
characteristics.  Like  all  magnetostrictive  materials,  its  strain 
(<)  is  a  function  of  both  the  stress  (a)  and  the  applied 
magnetizing  field  (H)  as 

«  -  f(o)  +  9(H)  (32) 

This  relation  for  Terfenol-D  is  shown  in  Figure  28.  Shown  are 
curves  of  strain  (in  parts  per  million  or  microstrain)  versus 
applied  magnetizing  field  (in  Oersteds)  for  various  levels  of 
compressive  stress.  These  curves  are  found  placing  a  Terfenol-D 


Applied  Field  (Oersted) 


Figure  28  Magnetostriction  of  Terfenol-D  vs 
Applied  Magnetic  Field 
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rod  inside  an  electrical  solenoid  and  axially  compressing  the 
Terfenol-D  with  a  fixed  load  (fixed  stress).  The  electric 
current  in  the  solenoid  is  then  varied,  which  varies  the  applied 
magnetizing  field  (H) .  The  resulting  axial  expansion  (tensile 
strain)  of  the  Terfenol-D  rod  is  then  measured. 

Many  of  the  important  characteristics  of  Terfenol-D  can  be 
seen  in  Figure  2823.  The  maximum  strain  is  about  0.2%  (2000 
microstrain) ,  approximately  10  times  better  than  piezoelectric 
materials.  The  strain  curves  are  symmetrical  about  the  zero 
applied  field  point  reflecting  the  fact  that  Terfenol-D  contracts 
for  both  positive  and  negative  magnetic  field  orientations.  The 
magnetostrictive  phenomenon  is  much  stronger  when  the  material  is 
under  compressive  stress.  The  change  in  strain  for  unit  change 
in  applied  field  is  strongest  at  applied  fields  of  approximately 
500  Oersted.  Because  of  these  properties,  Terfenol-D  is  best 
utilized  when  its  is  preloaded  with  compressive  stress  and  biased 
with  a  dc  magnetic  field. 

The  operation  of  Terfenol-D  actuators  can  be  explained  with 
the  simple  subsystem  shown  in  Figure  29.  Shown  is  a  Terfenol-D 
rod  surrounded  by  an  electric  solenoid.  A  simple  model  of  this 
subsystem  would  Include  the  solenoid  coil  current  and  stress  in 
the  Terfenol-D  as  external  inputs.  The  outputs  are  the  strain  in 
the  Terfenol-D  rod  and  the  coil  voltage.  In  addition  to  this 
Terfenol-D  subsystem,  a  complete  model  of  the  actuator  would 
include  the  stress-strain  relation  that  the  mechanical  subsystem 
of  the  actuator  would  impose  on  the  Terfenol-D. 
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Terfenol-D  Rod 


—  Circumferentially 
Wound  Solenoid 

T 


Figure  29  Simple  Maqnetostrictive  Actuator 

The  mechanical  power  produced  by  the  actuator  can  be  found 
from  the  stress-strain  curve  over  a  cycle.  For  example,  if  the 
stress  and  applied  field  waveforms  of  Figures  30  and  31  are 
applied  to  the  Terfenol-D,  the  resulting  stress-strain  relation 
shown  in  Figure  *2  can  be  found  using  the  curves  of  Figure  28. 
The  area  enclosed  by  the  stress-strain  loop  of  Figure  32  is 
simply  the  mechanical  energy  per  unit  volume  that  is  produced  by 
the  Terfenol-D.  This  energy  density  can  be  as  high  as  25kJ/m3 
(0.67  BTU/ft3),  approximately  an  order  of  magnitude  higher  than 
for  piezoelectric  materials. 
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Applied  Magnetic  Field  (Oersted)  Stress  (MPa) 
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Figure  30  Example  Stress  Waveform  in  Terfenol-D 
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Figure  31  Example  Applied  Magnetic  Field  Waveform  in  Terfenol-D 
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Figure  32  Stress  vs  Strain  in  the  Terfenol-D  for  the 
Example  gave forma 

The  mechanical  power  output  of  the  actuator  depends  both  on 
the  energy  delivered  during  the  cycle  and  the  frequency  of 
operation.  If  the  actuator  were  run  at  1000  Hz,  power  densities 
of  up  to  2.5  kW/kg  (1.5  hp/lba)  are  possible.  Note  that  this  is 
approximately  10  times  the  power  density  (approximately  o.l 
hp/lbjn)  for  conventional  electric  motors.  Higher  frequency 
operation  resulting  in  higher  power  densities  is  also  possible. 
The  efficiency  (mechanical  output  power/electrical  input  power) 
of  the  actuator  can  be  50%  or  greater.  Losses  are  dominated  by 
the  magnetic  hysteresis  losses  in  the  Terfenol-D,  but  also 
include  resistive  losses  in  the  solenoid  and  eddy  current  losses, 
which  can  be  minimized  by  the  use  of  laminations,  in  the 
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Terfenol-D. 

Besides  the  stress/strain/magnetic  field  properties  shown  in 
Figure  28,  other  important  properties  of  Terfenol-D  are  given  in 
Table  3  below25.  The  high  compressive  strength  of  Terfenol-D 
combined  with  the  greater  magnetostriction  when  under  compression 
means  that  Terfenol-D  should  nominally  be  under  compression. 


Table  3.  Terfenol-D  Properties 


Msshanlcfll 

Density 

Young's  Modulus 
Tensile  Strength 
Compressive  Strength 
Thermal 

Thermal  Expansion 
Magnetic 


9.25  103  kg/m3 
25-35  GPa 
28  MPa 
700  MPa 

12  X  10 


0.33  lb*/ in3 
3. 6-5.0  10®  psi 
4.1  ksi 
loo  ksi 


6  /degree  C 


Relative  Permeability  4-6 

Mechanical  Coupling  Factor  0.7  -  0.75 


In  summary,  Terfenol-D  is  an  exciting  new  material  that  will 
allow  the  design  of  novel  electromagnetic  actuators.  It  is 
particularly  well  suited  to  applications,  such  as  distributed 
actuators,  that  require  a  combination  of  large  forces,  relatively 
small  motions,  high  bandwidth,  and  good  transient  response. 
Compared  to  piezoelectric  materials,  it  offers  a  factor  of  ten 
improvement  in  strain,  energy  density,  and  power  density  without 
requiring  the  high  voltages  typical  of  piezoelectric  actuators. 
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Properly  designed,  Terfenol-D  based  actuators  can  hava  power 
dansitiss  an  order  of  magnitude  greater  than  conventional 
electric  motors  and  solenoids.  In  addition  Terfenol-D  based 
actuators  can  be  simply  controlled. 

This  brief  investigation  into  the  properties  of  Terfenol-D 
show  that  exciting  new  actuator  designs  can  be  developed  based  on 
its  unique  properties. 

Summary  of  Actuator  Design  Section 

In  this  section,  magnetic  analyses  were  reviewed  briefly  in 
order  to  present  the  actuator  design  types  considered.  Five 
actuator  designs  were  considered:  interaction  of  distributed 
coils  with  the  Earth's  magnetic  field,  interaction  of  distributed 
coils  with  a  large  flux  source,  interaction  of  distributed  coils 
with  each  other ,  local  interaction  of  coils  on  a  flexible 
structure  truss  element,  and  use  of  a  magnetostrictive  material, 
Terfenol-D.  The  first  two  approaches  were  shown  to  be 
undesirable  for  control  of  the  beam  model  used  in  this 
investigation.  Further  comparative  analysis  of  the  distributed 
coil  interaction  approach  and  the  local  interaction  approach  is 
presented  in  the  next  section.  The  use  of  Terfenol-D  is  very 
promising  and  requires  a  more  complete  analysis.  This  analysis 
is  being  proposed  as  an  initial  task  for  Phase  II. 

COMPARISON  OF  ACTUATOR  APPROACHES 
In  this  section  further  analysis  is  presented  to  compare  the 
more  promising  actuator  approaches,  ie  the  distributed  coil 
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interaction  and  the  local  coil  interaction.  The  actuator 
approach  requiring  the  use  of  Terfenol-D,  which  is  very 
promising,  requires  further  analysis  in  Phase  XX.  The  local  and 
distributed  actuator  approaches  were  coapared  in  two  ways.  The 
first  comparison  was  a  mass  and  power  comparison  and  was 
calculated  by  Mtching  actuator  force.  The  second  comparison  was 
done  on  the  basis  of  energy  removed  from  the  baseline  flexible 
structure  model  developed  under  task  1.  The  second  methodology 
quantifies  both  the  concepts  of  control  effectiveness  and  the 
optimum  actuator  location  for  control  of  a  vibrating  beam.  The 
control  effectiveness  is  a  function  of  both  actuator  type  and 
mode  shape . 

Force  Comparison 

One  way  to  compare  the  local  and  distributed  actuator 
schemes  is  to  match  the  vertical  force  transmitted  by  the 
actuators  to  the  beam  and  then  to  compare  the  required  mass  and 
power  to  obtain  this  force  (For  the  local  interaction  scheme,  the 
truss-beam  interaction  is  modelled  as  if  the  beam  imposed  a 
kinematic  constraint  of  no  deflection,  thus  force  and  'not 
deflection  is  transmitted  to  the  beam.).  For  the  local  actuator 
approach,  mass  and  power  comparison  calculations  were  performed 
using  a  current  density  of  10*  A/m2  (9  x  106  A/ft2) ,  wire  cross- 
sectional  area  of  1  cm2  (0.155  in2),  and  coil  radius  of  1.27  cm 
(1/2  inch) .  The  ratio  of  separation  distance  (r)  to  coil  radius 
(#  ■  1.27  cn  ■  1/2  in)  is  varied  from  1  to  4.  This  is  the  near 
field  range  where  Maxwell's  mutual  inductance  formula  (30)  is 
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appropriate.  Those  values  are  also  appropriate  tor 
implementation  in  the  HAST  model.  Aluminum  conductor  wee  assumed 
for  both  the  local  and  the  distributed  actuator  approaches, 
parameters  compatible  with  the  MAST  model  were  also  chosen  for 
the  distributed  case.  In  this  case,  coil  separation  distance  (r) 
was  equal  to  the  distance  between  the  distributed  coils  if  they 
were  located  at  the  ends  of  a  MAST  truss  section.  Because  the 
MAST  trusses  were  1.6  m  long  (5.2  ft)  with  a  truss  angle  of  47°, 
the  separation  distance  for  the  coils  is  1.6cos47°  m  (5.2  cos47° 
ft) .  The  coil  radius  for  the  coils  in  the  distributed  case  were 
chosen  to  be  5  cm  (1.97  in).  Notice  that  this  value  of  r/a  *  20 
is  appropriate  for  the  dipole  far  field  equations  (Equations  (25) 
and  (26)).  Magnetic  moments  of  the  coils  for  the  distributed 
case,  were  chosen  to  match  the  vertical  force  transmitted  to  the 
beam  calculated  in  the  local  case.  Figure  33  is  a  plot  of  mass 
vs  vertical  force  transmitted  to  the  beam.  Figure  34  is  a  plot 
of  power  vs.  vertical  force  transmitted  to  the  beam.  If  the 
coils  were  superconducting,  no  power  would  be  required.  The  plots 
show  that  the  local  interaction  scheme  requires  both  less  mass 
and  less  power  than  the  distributed  scheme.  The  larger  forces 
are  obtained  for  the  cases  where  the  coils  are  close  together. 
Thus  r/a-1  is  the  point  corresponding  to  the  largest  force  of 
approximately  0.4  N  (0.09  lb)  in  Figures  33  and  34.  This 
comparison  is  useful  as  it  gives  the  tradeoff  between  mass, 
power,  and  force. 

*• 
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Muss  Comparison  for  Local  vs.  Distributed  Interaction  Schemes 


Vertical  Force  on  Beam  (Newtons) 


x  Local 
o  Distributed 


Figure  33 


Control  Effectiveness  and  Optima  Actuator  Location  Deriva¬ 
tion  of  Control  Effectiveness 

A  methodology  is  developed  in  this  section  to  determine  the 
optimum  location  for  removing  power  from  a  vibrating  beam  using 
the  local  or  distributed  actuators.  The  local  and  the 
distributed  coil  interaction  schemes  are  then  compared  using  this 
methodology.  The  approach  is  to  consider  the  energy  removed  from 
the  beam  as  a  function  of  actuator  force  and/or  torque  and  modal 
displacements,  and  then  maximize  the  energy  removed  with  respect 
to  the  modal  parameters. 
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x  Local 
o  Distributed 

Figure  34  Power  Comparison  of  Local  vs  Distributed 
Actuator  Concepts 


For  example,  an  actuator  which  produces  a  single  force  would 


remove  power  from  the  beam  as  follows: 

P  *  F*y»« 


where 


(32) 


F  -  Actuator  Force 


y  =  Translational  Displacement  at  Actuator  Location 
u  ■  Modal  Frequency 

Therefore  the  energy  removed  per  cycle  is  given  by: 

E  -  F*y  (33) 

On  the  other  hand,  an  actuator  which  produces  a  single  torque 
would  remove  power  from  the  beam  given  by  the  following 
expression.  " 

P  =  r*  e.  w  (34) 

where 

r  =  Actuator  Torque  ‘ 

e  =  Rotational  Displacement  at  Actuator  Location 
w  *  Modal  Frequency 

Energy  removed  per  cycle  by  a  single  torque  actuator  is 
therefore: 

E  *  r • 8  (35) 

The  actuators  considered  here  are  not  single  force  or  torque 
actuators.  They  produce  both  force  and  torque.  A  general 
expression  can  be  written  for  the  energy  removed  per  cycle  from  a 
beam: 

n  m 

E  =  Si  Fi  yi  +  Sj  rj  9j  (36) 

1  1 

where  there  are  n  actuator  forces  and  m  actuator  torques.  This  is 
obtained  by  summing  the  energy  removed  per  cycle  at  each  actuator 
location.  This  term  will  be  referred  to  in  this  paper  as  the 
control  effectiveness  of  the  actuators.  Later  comparisons  of  the 
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ability  of  the  different  actuator  types  in  damping  out  the 
vibrations  of  a  beam  will  be  made  using  this  concept. 

If  a  sinusoidal  mode  shape  is  assumed,  the  translational 
displacement  (y)  for  the  beam  can  be  written  as: 

2  «  x 

y  -  A  sin  ( - )  (37) 

1 

where  x  is  the  axial  position  and  1  is  the  beam  length. 

Using  Bernoulli-Euler  beam  theory  the  rotational  displacement, 
e,is  obtained  by  taking  the  derivative  with  respect  to  x. 


2  it  2  f  X 

e  «  (  —  )  cos  ( - )  (38) 

1  1 


Comparison  of  Control  Effectiveness  for  Distributed  and  Local 
Cases 

Expressions  for  the  energy  removed  per  cycle  (effectiveness) 
will  be  derived  for  the  distributed  coil  interaction  case  and 
several  variations  of  the  local  coil  interaction  case. 

The  truss  beam- interaction  for  the  local  case  is  modelled  as 
if  the  beam  imposes  a  kinematic  constraint  of  no  beam  deflection. 
This  type  of  model  was  shown  in  Figure  27(a)  where  the  force 
transmitted  from  the  truss  to  the  beam  is  shown  to  be 
proportional  to  the  ratio  of  coil  separation  distance  to  truss 
length.  There  are  three  variations  of  the  local  actuator 
approach  that  are  considered.  The  first  is  a  one  truss  actuator 
that  is  in  tension  (or  compression) .  The  second  is  a  two  truss 
actuator  in  which  both  trusses  are  simultaneously  either  in 
tension  or  in  compression.  The  third  approach  is  a  two  truss 
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actuator  in  which  one  truss  is  in  tension  and  the  other  truss  is 
in  compression.  These  variations  on  the  local  actuator  approach 
and  associated  force  modelling  assumptions  are  shown  in  Figure 
35.  Notice  from  these  figures  that  axial  effects  on  the  beam  are 
neglected  because  only  beam  bending  is  being  modelled.  These 
models  are  appropriate  for  use  in  conjunction  with  the  finite 
element  beam  model  developed  previously.  The  elements  are  point 
nodes;  hence,  the  model  of  the  beam-truss  interaction  is  one 
which  uses  an  equivalent  force  system  on  a  node  located  on  the 
neutral  axis  of  the  beam.  Notice  in  Figure  35  the  notation  ••Ft" 
denotes  the  force  transmitted  to  the  beam  by  the  truss,  ”e"  is 
the  truss  angle,  and  ”un  is  the  beam  width.  The  one  truss 
actuator  produces  a  moment  and  a  force  on  nodes  which  correspond 
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(a)  One  Truss  Actuator 

Modelling  Assumptions  Appropriate  for  Use  in 
Finite  Element  Beam  Model  for  Local  Interaction 
Actuator  Schemes  (continued  on  next  page) 
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(b)  Two  Truss  Actuator  in  which  Both  Trusses 
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(c)  Two  Truss  Actuator  in  which  One  Truss  is  in 
Tension  and  the  other  is  in  Compression 


Figure  35 


to  the  axial  location  of  the  ends  of  the  truss  Figure  35(a).  The 
force  and  torque  on  the  structure  produced  by  this  actuator,  as 
well  as  the  other  types  of  actuators,  is  zero,  as  expected  from 
Newton's  Action-Reaction  Law.  The  two  truss  actuators  can  be 
modelled  by  force  and  torque  effects  at  three  nodes.  For  the  two 
truss  actuator  that  has  both  trusses  in  either  tension  or 
compression  (Figure  35(b)),  the  two  outer  nodes  undergo  forces  of 
the  same  direction  and  magnitude  and  torques  of  opposite 
directions  and  the  same  magnitude.  The  middle  node,  which 
corresponds  to  the  truss  intersection  location,  undergoes  a  force 
that  is  twice  the  magnitude  of  one  of  the  individual  forces  on  an 
outer  node.  Again,  by  adding  up  these  forces  and  torques,  the 
net  effect  is  zero.  The  third  actuator  approach  (two  truss 
compression-tension)  produces  a  different  effect  on  the  overall 
beam  structure  (Figure  35(c)).  The  middle  node  has  a  moment  that 
is  twice  as  strong  and  in  the  opposite  direction  as  the  moments 
on  the  outer  nodes.  The  outer  nodes  have  forces  which  are  equal 
in  strength  but  opposite  in  direction. 

The  distributed  actuator  approach  can  be  modelled  in  a 
similar  way  (Figure  19) .  This  model  has  two  nodes  located  at  the 
positions  of  the  two  coils.  The  two  forces  produced  are  of  equal 
magnitude  and  opposite  direction.  The  two  torques  produced  in 
the  coils  cancel  the  net  moment  produced  by  the  forces. 

Because  the  distributed  case  and  the  one  truss  local  case 
reduce  to  similar  force  effects  on  the  beam,  a  common  terminology 
is  introduced  (Figure  36) .  The  mean  location  of  the  actuator 
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by  the  wavelength  of  the  mode. 

n  1 

P  »  - 

*b 

where  1  Is  the  distance  between  force  application  points. 
Expressions  for  the  control  effectiveness  E  (energy  removed  per 
unit  cycle)  were  derived  using  Equation  (36)  for  each  actuator 
type  being  compared.  These  expressions  are  shown  in  Table  4. 


Tabic  4  Effectiveness  Expressions  for  Local  and 
Distributed  Actuator  Types 


Local  One  Truss 


Effectiveness  Expression 
2FA  cos  a [0  cos  p  -  sin  p ] 


Local  Two  Truss 
Tension-Tension 


_  p 

4FA  sin  a[sin20  -  -  sin  (20)] 
small  p  E  a  04 


Local  Two  Truss  2AF  cos  a [sin  (29)  -  2 p  cos2  0] 

Compression-Tension 

small  0  E  a  03 


Distributed  "Configuration 
2"  Scheme 


2FA  COS  a  sin  0  -  20FA  cos  a  cos  0 
2 

+  ^  0FA  sin  a  sin  0 


In  these  expressions,  "A"  is  the  amplitude  of  vibration  and 
"F"  is  the  vertical  force  transmitted  to  the  beam.  Effectiveness 
was  calculated  in  order  to  determine  which  local  actuation  scheme 
is  the  most  useful  for  flexible  structure  control  of  the  MAST 
beam  model.  The  calculations  were  performed  for  one  actuator 
located  at  its  optimum  location  on  the  flexible  beam  model 
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developed  previously.  The  opt  1  nun  location  Is  found  by 
maximising  the  effectiveness  term  with  respect  to  the  mean 
location  naN.  The  results  of  these  calculations  are  shown  in 
Figures  37  and  38.  Figure  38  is  plotted  on  a  logarithmic  scale 
to  emphasise  what  happens  when  is  small.  This  is  of  interest 
because  the  flexible  structures  have  low  natural  frequencies, 
long  wavelengths  and  thus  small  0's.  Using  Figure  38,  the  local 


Model  Control  Effectiveness  vs.  Coil  Separation 


Coil  Separation  (fraction  of  modal  wave  length) 


CASE  lx  1  TRUSS  EXTENSION 

CASE  2:  2  TRUSS  EXTENSION-EXTENSION 


CASE  3:  2  TRUSS  COMPRESSION-EXTENSION 


Normalized  Modal  Control  Effectiveness 


interaction  variations  can  be  compared.  This  plot  ahovs  that  an 
actuator  consisting  of  two  trusses,  one  undergoing  compression 
and  the  other  undergoing  extension,  is  Jiore  effective  for  small 
p  'a  than  a  one  truss  extension  or  a  two  truss  extension- 
extension.  In  fact,  it  can  be  shovn,  using  Taylor  series 
expansions,  that  the  effectiveness  for  the  truss  extension- 
compression  case  is  proportional  to  while  the  effectiveness 
for  the  truss  extension-extension  case  is  proportional  to  /94. 

Table  5  is  a  summary  of  effectiveness  calculations  on  a  per 
mass  and  per  power  basis  for  the  distributed  interaction  scheme 
and  two  variations  on  the  local  one  truss  interaction  scheme. 
Because  the  MAST  deployable  beam  structure  has  truss  angles  of 
47°  for  some  of  its  structural  trusses,  one  local  scheme 
considered  in  these  comparisons  also  uses  a  truss  angle  of  47°. 
Analysis,  however,  showed  that  for  the  purpose  of  removing 
energy  from  a  vibrating  beam  with  the  actuators  considered,  a 
truss  angle  of  30°  is  more  effective.  Thus  a  second  local  scheme 
is  considered  with  a  truss  angle  of  30°.  The  modal  wavelength 
used  in  these  calculations  is  60  m  (197  ft) ,  which  is  the  length 
of  the  MAST  structure,  and  a  non-superconducting  current  density 
level  of  106  A/m2  (9  x  104  A/ft2)  is  assumed.  The  chart  shows 
that  these  local  interaction  schemes  are  two  orders  of  magnitude 
better  than  the  distributed  interaction  scheme  on  both  an 
effectiveness  per  mass  and  an  effectiveness  per  power  basis.  The 
best  effectiveness  per  mass  per  vibration  amplitude  (energy 
removed  per  cycle  from  a  vibrating  beam  per  Kg  of  actua'cor  mass 
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Table  5  Effectiveness  Calculations 
Modal  Wavelength:  60  a 
Current  Density:  10®  A/®2 

E  J  E  J  E  jJ 

A  a  cycle  An  nXg  cycle  AP  n  Watt  cycle 


Distributed  7.78X10”10  2.38  x  10“10  2.6  x  10“H 

Interaction 

"Standard"  Sise 

Dipoles  10  m  Apart 

One  Truss  Local  3.2xl0“10  2.9  x  10“8  2.9  x  10“9 

Interaction  Scheme 

for  MAST 

Truss  Angle  47° 

Truss  Length  1.6  m 
Wire  C.S.  Area  (l/2)2cro2 
Coil  Radius  1.27  cm 

One  Truss  Local  4.5x10-1°  4.1  x  10“8  4.1  x  10“9 

Interaction  Scheme 
Truss  Angle  30° 

Truss  Length  1.6  m 
Wire  C.S.  Area  (1/2) 2cm2 
Coil  Radius  1.27  cm 


per  vibration  amplitude)  in  this  chart  is  4  x  10~8  J/m  kg  cycle 
(5  x  10-12  etu/ft  lbm  cycle)  and  4  xlO”9  J/m  watt  cycle  (8.6  x 
10-10  Btu/ft  hp  cycle)  for  the  local  two  30°  truss  extension- 
compression  case. 

Because  of  the  simplifying  modelling  assumptions  used  in 
this  Phase  I  investigation,  these  effectiveness  numbers  should 
only  be  used  for  comparing  the  different  actuator  approaches.  In 
particular,  these  calculations  are  based  on  extreme  models  of 
truss-beam  interactions.  In  either  extreme  model  discussed 
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earlier  (aero  foroe  va  saro  deflection  Figure  27),  tha  desired 
affact  (forca  or  daflaction)  ia  raducad  by  a  factor  of  r/1 
(aaparation  diatanca  to  truaa  langth)  if  tha  coila  ara  placed  in 
tha  middle  of  tha  truaa.  If  two  aata  of  coila  ara  placad  cloae 
to  tha  anda  of  tha  truss,  however,  an  equivalent  simple  model 
will  increase  by  a  factor  of  1/r  tha  effectiveness  calculated 
using  tha  rigid  beam  modal. 

Tha  coils  are  capable  of  producing  adequate  forces.  As 
shown  earlier,  with  current  densities  of  108  A/m2  (9  x  106 
A/ ft2) ,  the  local  interaction  case  was  shown  to  be  capable  of  5  N 
(1.1  lb)  for  geometrical  parameters  appropriate  for  MAST. 
Because  of  the  simplifying  modelling  effects,  further 
effectiveness  analysis  is  suggested  as  an  initial  task  in  Phase 
II. 

Tables  6  and  7  are  similar  charts  which  show  the  same  effect 
for  different  wavelengths  and/or  current  densities.  The  trend 
showing  the  local  interaction  is  better  than  the  distributed 
interaction  in  an  effectiveness  per  mass  basis  is  again  shown  in 
these  charts.  Table  6  is  for  a  wavelength  of  360  m  (1181  ft), 
and  a  non-superconducting  current  density.  By  comparing  this 
chart  with  the  previous  chart  (Table  5)  it  can  be  seen  that  the 
effectiveness  decreases  with  increasing  wavelength.  This  is  to 
be  expected  as  the  equations  for  effectiveness  for  small  p  (Table 
4)  show  that  the  effectiveness  is  proportional  to  a  power  of  p 
and  thus  inversely  proportional  to  wavelength. 
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Table  6  Effectiveness  calculations 
Modal  Wavelength:  360  a 
Currant  Oanaity:  10®  A/a2 

E  J  E  J  E  J 

A  a  cycla  Aa  aXg  cycla  AP  a  Matt  cycle 


Distributed  2.0X10”11  6.1  x  10“12  6.1  x  10“13 

Interaction 

"Standard"  Size 

Dipoles  10  n  Apart 

One  Truss  Local  1.5xl0“12  1.3  x  10“30  1.3  x  lO-*11 

Interaction  Scheme 

for  MAST 

Truss  Angle  47° 

Truss  Length  1.6  m 
Wire  C.s.  Area  (l/2)2cm2 
Coil  Radius  1.27  cm 

One  Truss  Local  2.lxl0“12  1.9  x  10“10  1.8  x  10*"11 

Interaction  Scheme 
Truss  Angle  30° 

Truss  Length  1.6  m 
Wire  C.S.  Area  (1/2) 2cm2 
Coil  Radius  1.27  cm 


Table  7  is  effectiveness  calculations  using  a 
superconducting  current  density  wavelength  of  360  m  (1181  ft) , 
respectively.  The  effectiveness  per  power  basis  is  not  shown 
here  because  superconductors  use  no  power.  These  numbers  are  an 
improvement  by  4  orders  of  magnitude  over  their  respective 
effectivenesses  in  Tables  6.  The  actuator  type  with  the  highest 
effectiveness  can  remove  1.9  x  10”6  J/m  fcg  cycle  (2.5  x  10“10 
Btu/ft  lbm  cycle) . 
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Table  7  Effectiveness  Calculations 


Modal  Wavelength:  360  n 
Current  Density:  10®  A/m2 


E  _ 

A  m  cycle 


Distributed  2.0  x  10"7 

Interaction 
"Standard”  Size 
r Ipoles  10  m  Apart 

One  Truss  Local  1.5  x  10“8 

Interaction  Scheme 

for  MAST 

Truss  Angle  47° 

Truss  Length  1.6  m 
Wire  C.S.  Area  (1/2) 2 cm2 
Coil  Radius  1.27  cm 

One  Truss  Local  2.1  x  10"® 

Interaction  Scheme 
Truss  Angle  30° 

Truss  Length  1.6  m 
Wire  C.S.  Area  (1/2) 2cm2 
Coil  Radius  1.27  cm 


E_  __  j__ 
Am  mKg  cycle 


6.3  X  10“8 


1.3  X  1CT6 


1.9  X  10-6 


Summary  of  Actuator  Comparison  Section 

In  this  section  two  types  of  actuator  designs  were  compared: 
the  local  coil  interaction  and  distributed  coil  interaction 
cases.  Comparisons  were  performed  on  both  a  force  and  an  energy 
basis.  The  local  coil  interaction  was  shown  to  be  a  better 
approach  on  a  per  mass  and  a  per  power  basis  for  controlling  beam 
vibrations.  In  addition,  several  variations  on  the  local  coil 
interaction  scheme  were  examined. 
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CONTROLLER  DESIGN 


The  purpose  of  this  section  is  to  verify  the  open  loop 
analysis  developed  previously  for  the  effectiveness  comparison  of 
the  various  actuator  schemes.  This  will  be  done  by  simulating 
the  MAST  model  with  approximately  200  local  actuators  of  five 
Newton  force  (1.1  lb)  capability  actuators  placed  optimally  on 
the  structure.  As  stated  earlier,  the  effectiveness  numbers  are 
strongly  dependent  on  the  modelling  assumptions  of  the  truss-beam 
interactions  and  should  only  be  used  to  compare  the  different 
actuator  approaches. 

The  method  used  to  verify  the  open  loop  predictions  is  to 
plot  time  simulations  of  a  vibrating  beam  placed  under  optimal 
control  action  using  two  different  local  actuation  schemes.  The 
time  simulations  should  show  that  the  scheme  with  the  higher 
effectiveness  will  be  able  to  damp  the  vibrations  faster. 

Two  local  interaction  schemes  are  shown  along  with  a  diagram 
of  the  forces  and  moments  which  are  transmitted  to  the  beam. 
Each  local  interaction  scheme  consists  of  two  truss  structures. 
These  are  shown  in  Figures  35(b)  and  35(c).  Each  truss  has  a  pair 
of  coils  centered  in  the  truss  section.  In  the  first  local 
interaction  scheme  the  coil  forces  are  set  up  to  produce  tension 
in  both  trusses  or  compression  in  both  trusses.  The  second  local 
interaction  scheme  is  set  up  to  produce  compression  in  one  truss 
and  tension  in  the  other  truss.  Time  simulations  were 
performed  with  a  mode  3  initial  condition  and  a  maximum 
transverse  displacement  of  1  cm  (0.4  in).  This  mode  has  a 
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wavelength  that  is  approximately  the  length  of  the  beam  (60  m 
(197  ft))  and  a  frequency  of  5.7  Hz.  The  maximum  transverse 
displacement  occurs  at  the  20  m  (65  ft)  node.  The  mode  shape  is 
shown  in  Figure  5.  Time  simulations  of  the  controlled  structure 
are  shown  in  Figures  39  and  40.  The  maximum  vertical  force 
transmitted  to  the  beam  was  matched  for  both  cases  to  be  around 
1000  N  (225  lb) .  The  effectiveness  calculations  presented 
previously  predict  that  the  compression  extension  scheme  will 
remove  more  energy  per  unit  time  from  a  vibration  beam  than  the 
extension-extension  case.  The  time  simulations  verify  this 
prediction.  They  show  higher  damping  levels  for  the  compression 
extension  scheme. 

CONCLUSIONS 

This  Phase  I  SBIR  research  has  investigated  the 
feasibility  of  using  an  array  of  magnetic  coils  for  use  as  an 
actuation  system  for  flexible  structure  control.  Five  actuator 
designs  iwere  compared.  The  first  actuator  design  involved 

t 

magnetic  coils  distributed  in  a  flexible  structure.  These  coils 
would  produce  control  torques  on  the  beam  by  interacting  with  the 
Earth's  magnetic  field.  The  second  actuator  design  was  control 
torque  production  by  the  interaction  of  distributed  magnetic 
coils  with  a  large  flux  source  located  on  the  spacecraft.  The 
third  actuator  design  produced  control  forces  and  torques  by 
interacting  distributed  coils  on  the  structure.  The  fourth  was 
to  place  the  coils  very  close  together  to  obtain  stronger  forces 
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and  torquas  on  the  flexible  structure.  This  approach  was  termed 
local  interaction.  The  final  actuator  approach  requires  the  use 
of  a  new  material,  Terfenol-D,  which  is  maqnetostictive ,  i.e.,  it 
strains  under  the  presence  of  a  magnetic  field. 

Initial  investigations  on  each  of  these  methods  wrre 
performed.  The  interaction  of  distributed  coils  with  the  Earth's 
magnetic  field  would  require  the  shape  and  vibration  control 
system  to  depend  strongly  on  the  position  of  the  actuator  with 
respect  to  inertial  space.  Also,  this  method  can  only  produce 
control  torques  because  the  force  depends  on  the  gradient  of  the 
magnetic  field  which  is  small  for  the  distances  associated  with 
the  magnetic  coils.  The  second  method,  interaction  with  a  large 
flux  source,  is  only  capable  of  torque  production  about  two  axes. 
And  furthermore,  there  would  be  a  torque  on  the  large  flux  source 
which  is  large  enough  to  interfere  with  attitude  control.  These 
first  two  actuator  methods  are  undesirable  for  a  HAST  type 
structure.  They  may  have  application,  however,  to  a  spacecraft 
with  combined  rigid  body  attitude  control  and  flexible  structure 
control.  The  actuator  approach  using  a  new  magnetostrictive 
material,  Terfenol-D,  was  shown  to  be  a  very  promising  approach 
that  requires  further  Phase  II  investigation. 

The  emphasis  of  this  Phase  I  research  was  comparing  the 
distributed  coil  and  the  local  coil  interaction  methods.  The 
local  interaction  approach  transmits  more  force  to  the  structure 
on  both  a  per  mass  and  a  per  power  basis  than  does  the 
distributed  coil  approach.  A  concept  was  developed  which 
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quantifies  tha  effectiveness  of  an  actuator  being  used  to  damp 
vibrations  in  tha  MAST  beam  modal.  This  concept  is  based  on  the 
amount  of  energy  removed  per  cycle.  Several  variations  of  the 
local  interaction  actuator  concept  were  then  investigated.  These 
included  a  one  truss  local  actuator  and  two  truss  actuators.  One 
type  of  two  truss  actuator  had  tension  (or  compression)  in  both 
trusses.  The  other  type  had  tension  in  one  of  the  trusses  and 
compression  in  the  other  truss. 

The  open-loop  effectiveness  comparisons  were  verified  by 
simulation.  A  controller  was  designed  for  the  system  using 
linear-quadratic  methods.  The  two  truss  tens ion-tens ion  case  and 
the  two  truss  compression-tension  case  were  compared  and  verified 
the  control  effectiveness  analysis. 

Further  comparative  analysis  of  the  local  interaction 
approach  and  the  Terfenol-D  approach  is  recommended  for  Phase  II 
research.  Simpler  analytical  models,  for  example  those  that 
model  strictly  the  actuator  and  truss  section,  are  recommended. 
Phase  II  will  include  the  design  and  construction  of  both  types 
of  prototype  actuators.  The  use  of  the  truss  and  actuator  models 
will  also  allow  comparison  with  more  conventional  actuator 
approaches  such  as  piezoelectrics. 
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